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We study the effect of heterogeneous load sharing in the fiber bundle models of fracture. We divide 
the system into two groups of fibers (fraction p and 1 — p) in which one group follow the completely 
local load sharing mechanism and the other follow global load sharing mechanism. We argue that 
patches of local disorders in the loading plate can cause such a situation. We find that in 2d a finite 
crossover point comes up near p c ~ 0.53, which is slightly below the site percolation threshold. We 
numerically determine the phase diagrams (in Id and 2d) and identify the critical behavior below 
p c with the mean field behavior (completely global load sharing). This crossover can occur due 
to geometrical percolation of disorders in the system and the critical point can get modified by 
occasional long range load transfer process. We also show how the critical point depends on the 
loading history, which is identified as a special property of local load sharing. 



I. INTRODUCTION 

The fiber bundle model serves as a simple model for fail- 
ure due to fracture. Originally starting with the tex- 
tile industry [l[ , some of its characteristic features mimic 
the failure processes in many disordered materials . 
Also, it can have direct similarities with fiber reinforced 
composites, which are very useful materials for aerospace 
industries [fj. Due to its importance in theoretical mod- 
elling of disordered materials in general and FRCs in par- 
ticular, fiber bundle model have been studied in great 
details 

The basic model is about a bunch of fibers, having dif- 
ferent failure thresholds and hanging from a rigid plate. 
The fibers are also attached to a plate at the bottom, 
from which a load is hanging. If the threshold distribu- 
tion starts continuously from zero (in most cases they 
do), then as soon as any load is applied, some of the 
weak fibers will break. The load carried by those fibers 
will now be redistributed among the other intact fibers 
and this may lead to failure of some more fibers and so on. 
This avalanche can stop due to one of the following two 
reasons: cither all the fibers may break and the system 
undergoes complete failure, or the remaining fibers can 
withstand the load and the system is only partially dam- 
aged. The transition between partial and complete fail- 
ure have been studied using the standard tools of phase 
transitions (diverging time scales, scale free behavior of 
avalanche size distribution, susceptibility and so on). 

Depending upon the requirements and interests of 
study, the fiber bundle model can differ in some key as- 
pects. The threshold distribution of the fibers can vary 
(uniform distribution, Gaussian, Wcibull, Gumbel etc.). 
Although the failure point sensitively depends upon the 
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details of these distributions, the exponents remain uni- 
versal. Another key aspect in which the models can differ 
is the load sharing mechanisms. When a fiber breaks, its 
load is to be shared by the remaining intact fibers. But 
there can be many ways in which this can be done. Two 
extreme cases are: Global Load Sharing (GLS) and Local 
Load Sharing (LLS). In GLS, the load of the broken fiber 
is distributed equally among the rest of the fibers and in 
LLS it is distributed only among the nearest surviving 
neighbor(s). In general, the elastic properties of the bot- 
tom plate (or that of the matrix material in FRCs) are 
responsible for the load sharing mechanisms. The two 
extreme cases correspond to perfectly rigid and perfectly 
soft bottom plates. The real situation is, of course, some- 
where between these two. Hence, there have been efforts 
to interpolate between the two mechanisms [l7l - l20j . Fi- 
nally, another aspect that we want to mention here is 
the different loading processes. There are several ways in 
which the load can be increased. One way is to apply a 
load and then wait for the avalanche to stop. Once the 
avalanche stops, increase the load upto the point where 
the weakest surviving fiber breaks. This will start an- 
other avalanche and so on. Another way is to increase 
the load by a constant amount every time an avalanche 
stops. And finally, one can also directly apply the re- 
quired load in the intact system. In GLS models, the 
loading mechanisms do not change the critical point. In 
some situations, this can be important, as we shall dis- 
cuss. 

In the different models, the load sharing mechanisms 
followed by all the fibers were same (global, local or some 
other intermediate rule). However, the governing factors 
of the load sharing mechanisms, like the elastic proper- 
ties of the bottom plate, can be modified due to local 
disorders, i.e., the elastic properties of the bottom plate 
may not be homogeneous. The heterogeneity introduced 
by local disorders can make the load sharing mechanisms 
heterogeneous as well. In this study we investigate the 
effect of such heterogeneity. We divide the fiber into two 
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groups, where one group follow GLS and the other follow 
LLS. The competing effects of these two mechanisms give 
rise to a cross-over behavior at a threshold value of the 
fraction. We determine the phase diagram and the criti- 
cal behavior numerically. Also, we show that the process 
of loading changes the critical point in this case. We ar- 
gue that spatial correlations in the presence of local load 
sharing gives rise to this change in critical point. 



II. HETEROGENEOUS LOAD SHARING 

As pointed out before, the load sharing process in the sys- 
tem is a manifestation of the elastic properties of the base 
plate. Hence, local disorders in that material may lead to 
heterogeneous load sharing processes. To see some effects 
of that we divide the system into two parts, a fraction p 
of the fibers follow completely local load sharing process 
and the rest 1—p fraction follow global load sharing pro- 
cess. In increasing the load we have followed the direct 
loading mechanism, i.e., we apply a load directly in the 
intact fiber. Comparisons with other loading mechanisms 
are discussed in the later part . 

We study this model both in Id and 2d (square lattice) 
geometry numerically. Since this is also an interpolation 
mechanism (p = corresponds to GLS and p = 1 corre- 
sponds to LLS), we compare the results with that of rcf 
[18j , where a fraction g of the load of the broken fiber is 
shared locally and the rest is distributed globally. Below 
we present the simulation results in one and two dimen- 
sions for this model. 



A. Results in one dimension 

In this case the fibers are arranged in a linear lattice. 
After applying a load on the intact fiber, we scan the 
lattice once and the fibers having threshold below the 
applied load are broken. This would create patches of 
broken fibers. Now, among them the fraction of load 
carried by the local load sharing fibers, are distributed to 
the nearest surviving neighbors and the rest is distributed 
globally among the remaining fibers. This may lead to 
further failures and so on. 

While we plot the fraction of surviving fibers with ex- 
ternal load, to accurately determine the critical point we 
also study the reduced fourth order Binder cumulant [2l| 

t?-i (c/4) m 

which has the property that its value is independent of 
system size at the critical point. This helps in determin- 
ing the critical point very accurately. In Fig. [1] we plot 
the fraction of surviving fibers with external load and for 
different system sizes (at p = 0.3 and p = 0.0). Then 
in Fig. [2] we plot the Binder cumulant for different sys- 
tem sizes (also for p — 0.3 and p = 0.0). The common 
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FIG. 1: The fraction of surviving fibers arc plotted with 
load per fiber for p = 0.3 both in Id (left curves; 
L = 10000, 25000, 50000, 100000) and 2d (right curves; 
L = 50, 100, 200, 400, 700, N = LxL). 
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FIG. 2: The Binder cumulants are plotted with load per 
fiber for p = 0.3 both in Id (left curves; 
L = 10000, 25000, 50000, 100000) and 2d (right curves; 
L = 50, 100, 200, 400, 700, N = LxL). 



crossing point, which signifies size independent value, is 
the critical point. In this way the phase boundary was 
determined. 

To determine the critical behavior, we plot the relax- 
ation time in Fig. [3] The relaxation time is expected to 
diverge near the critical point as r ~ (<r c — &)~ z , From 
the simulations we see that z = 0.50 ± 0.01, which is 
in agreement with GLS exact result z = 1/2. We have 
checked that these exponent values remain unchanged 
upto p = 0.95 and beyond that no power law variation is 
observed. Of course p = 0.95 is the crossover point in the 
phase boundary, beyond which the system fails as soon as 
any finite load is applied and there is no common cross- 
ing point of the Binder cumulant. This indicates that the 
mean field critical behavior is present upto the crossover 
point and beyond that it goes over to local load sharing 
behavior, where of course critical behavior is not present. 
We also note that the crossover point is almost upto the 
point of percolation. In one dimension, if the disorder 
has to percolate, its concentration will have to be unity. 
In this case that value is almost reached. 
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FIG. 3: Relaxation time (r) is plotted for one 
dimensional heterogeneous model with p = 0.3,0.5. A 
power- law variation of the form r ~ \o — o c \~ z was 
observed (for both side of the critical point), with 
z = 1/2, similar to the global load sharing model. 
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FIG. 4: The fraction of surviving fibers are plotted 
against load per fiber for (left to right) 
p = 0.0, 0.3, 0.5, 0.52, 0.6 and for 
L = 50, 100, 200, 400, 700 for each p in two dimensions. 
Where the transition point became sharper with system 
size for p < p Cl for p > p c the critical load goes to zero, 
signifying a crossover to LLS behavior. 



B. Results in two dimensions 

Similar to the case in one dimension, this heterogeneous 
load sharing mechanism can be followed in two dimen- 
sions as well. In this case, after each scan of the lattice, 
the clusters of the broken fibers are identified. Then the 
fraction of the load to be shared locally is distributed 
along the cluster boundary and the rest is redistributed 
globally. The fraction of surviving fibers with external 
load for various p values are shown in Fig. 0] As be- 
fore, one can determine the critical point by the common 
crossing point of the Binder cumulant for different sys- 
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FIG. 5: The Binder cumulants are plotted with load per 
fiber for p = 0.0, 0.3 and different system sizes 

(L = 50, 100, 200,400, 700) in two dimensions. The clear 
common crossing point gives a c in each case. 
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FIG. 6: The Binder cumulants are plotted with load per 
fiber for p = 0.55,0.60 and different system sizes in two 
dimensions. In sharp contrast with the case of p < p c , 
here we don't find any common crossing point for 
different sizes. 



tern sizes (see Fig. [5]). The critical point decreases with 
increasing fraction of the local load sharing fibers. But 
above a threshold value p c = 0.53 ± 0.01, there is no 
common crossing point of the Binder cumulant (see Fig. 
[6]). This suggests a crossover to local load sharing. This 
point is also slightly below the site percolation threshold 
(p c = 0.5927 . . . ) [22(. The critical behavior, determined 
from the avalanche size distribution (Fig. [7]) and diver- 
gence of relaxation time (Fig. [8]), is same as that of GLS 
model. The phase diagrams for one and two dimensions 
for this model is shown in Fig. |H1 



III. DEPENDENCE ON LOADING PROCESS 

As mentioned before, there can be several ways in which 
external load can be applied in the system. The two 
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FIG. 7: The avalanche size distribution is plotted for 
p = 0.0, 0.5, 0.7. It is seen that for p = 0.0, 0.5 the size 
distribution follow the same power law with exponent 

2.5 ± 0.01 and for p = 0.7 > p c , there is no scale free 
behavior. 
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FIG. 9: The phase diagrams for heterogeneous load 
sharing in one and two dimensions. The crossover points 
are at p c = 0.95 ± 0.01 and p c = 0.53 ± 0.01 respectively. 
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FIG. 8: Relaxation time (r) is plotted for two 
dimensional heterogeneous model with p = 0.3,0.5. A 
power- law variation of the form r ~ \cr — o c \~ z was 
observed (for both side of the critical point), with 
z = 1/2, similar to the global load sharing model. 



extreme cases are (i) increasing the load upto the point 
when the weakest fiber breaks and (ii) applying the de- 
sired load directly on the intact fiber. Between these two 
ways, one can also adapt a third possibility when the 
load is increased by equal amount in every step (after an 
avalanche stops). In GLS scheme, the fraction of surviv- 
ing fibers and for that matter the critical point would 
not depend on how one increases the load. However, 
the situation can change when there is local load sharing 
mechanism present (even partially). 

In the case of our model, we have studied the three 
ways of increasing the load and found that the critical 
point differs for finite p values. In Fig. [10] , we have 
plotted the fraction of surviving fibers for equal load in- 
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FIG. 10: The fraction of surviving fibers are plotted 
against external load for two loading processes: equal 
load increment with S = 0.0001 and direct loading, for 

p = 0.5 and 0.0. For p = 0.0, the two curves for two 
loading processes exactly fall on top of each other, while 
those for p = 0.5 are significantly different. 



crement and direct loading for p = 0.5. The results show 
that the critical point is significantly different for the 
two cases. We have also plotted the same cases for p = 
(GLS model) but it is clear that the load increment mech- 
anisms have no effect on the GLS version. This is proba- 
bly due to the fact that spatial correlations in avalanches 
are present in the local load sharing, which depends on 
the load increment method. But for GLS there is no 
spatial structure as such. 

In Fig. [TT]wc plot, for a given system size and p value, 
the fraction of surviving fibers for different loading meth- 
ods. We find that the weakest fiber breaking and the 
direct loading are two limits, and equal load increment 
method interpolates between the two. For very small 
value of the fixed increment rate (6 = 0.0001), the curve 
is almost identical with that in weakest fiber breaking. 
This is because, for 6 ~ 1 /AT, the equal load increment 
method should be effectively same as that of weakest fiber 
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FIG. 11: The fraction of surviving fibers are plotted 
against external load for different loading processes 
when p = 0.3. For equal load increment with S = 0.0001 
almost matches with the curve obtained for loading 
upto weakest fiber breaking. The equal step load 
increment with 5 = 0.0004, 0.001 are also shown. They 
are in between the S = 0.0001 curve and that due to 
direct loading. 



IV. COMPARISON WITH MIXED LOAD 
SHARING MODEL 

The present model is also an interpolation scheme be- 
tween LLS and GLS. In fact this is the heterogeneous 
version of the mixed load sharing model [lH , where after 
failure of every fiber, a fraction g of the load of a bro- 
ken fiber is distributed locally and the rest is distributed 
globally. Although the load sharing on average is same 
for the two model, the local fluctuations introduced in 
the present version, makes it significantly different from 
the earlier one. 

The relaxation time divergence of this model in two di- 
mensions (Fig. ri"2j) . shows the same critical behavior as 
that of the one dimensional version (and GLS model). In 
Fig. [T3]we compare the phase diagrams of the two models 
in one and two dimension, which shows the effect of het- 
erogeneity as the phase boundaries do not coincide. The 
phase boundary and critical behavior of MMLS model in 
one dimension was reported in ref. (l8| . Here we deter- 
mine those in two dimensions. 



V. SUMMARY AND DISCUSSIONS 
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FIG. 12: Relaxation time (r) is plotted for two 
dimensional homogeneous model [18| with p = 0.3, 0.5. 
A power-law variation of the form r ~ |c — ct c | _z was 
observed (for both side of the critical point), with 
z = 1/2, similar to the global load sharing model. 



breaking. 

Finally, we note that although the critical point is de- 
pendent on the load increment mechanisms for finite p 
values, the crossover point and the critical behavior be- 
low the crossover point is unaffected by the load sharing 
mechanisms. The results presented in the previous sec- 
tion are for direct loading, except for the avalanche size 
distribution, which is for weakest fiber breaking. 



There are two extreme cases of load sharing (GLS and 
LLS) in the fiber bundle models. Since load sharing 
mechanism is the only source of cooperatively in these 
models, in that sense the load sharing mechanisms deter- 
mine the spatial fluctuation in that cooperatively. Apart 
from the theoretical significances of the possible spatial 
correlations, the load sharing mechanism is of vital prac- 
tical importance, since it determines the strength of the 
material (irrespective of threshold distribution). 

While the completely global load sharing is idealised 
case, the completely local load sharing has no finite crit- 
ical point. Hence, there have been several interpolation 
mechanisms to match realistic situations [r7H20| . How- 
ever, in all such cases the load sharing methods is same 
throughout the sample. But here we argue that local 
disorders in the materials can make the load sharing spa- 
tially heterogeneous. In this regards, we have proposed 
a simple heterogeneous load sharing model that gives 
crossover behavior at a non-trivial concentration value 
of the disorder. 

We have studied this model in both one and two di- 
mensions. In one dimension, we have found that the 
critical disorder concentration at which the crossover oc- 
curs is p c = 0.95 ± 0.01, which is very close to the triv- 
ial percolation threshold in one dimension (unity). The 
phase diagram was determined using Binder cumulant. 
It is seen that the critical behavior, in terms of avalanche 
size distribution and relaxation time divergence, remains 
that of the GLS model upto the crossover point (see Fig. 
[3]). The same behavior is seen in two dimensions. How- 
ever, the crossover point is now p c = 0.53 ± 0.01, which 
is again close (and slightly below) the site percolation 
threshold in square lattice (0.5927 . . . ) [22j. The fact that 
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FIG. 13: The left panel shows the comparisons of the 

phase diagram of the heterogeneous and the 
homogeneous models in one dimension and the right 
panel shows the comparisons of the phase diagram of 
the heterogeneous and the homogeneous models in two 
dimension. 



the crossover occurs even before the percolation may be 
due to the fact that even for local load sharing, the load 
transfer may some times be long ranged due to the ab- 
sence of neighbors. 

We have also noted that a finite local load sharing frac- 
tion would lead to dependence on loading history. For a 
GLS model, the fraction of surviving fibers and the criti- 
cal point was independent of loading history. But as soon 



as a local load sharing fraction is introduced, these quan- 
tities become loading history dependent (see Fig. ITTT) . 
This is an important factor, because in a real material it 
is not possible to determine the microscopic load sharing 
mechanism directly. But this is a macroscopic manifesta- 
tion of the microscopic load sharing process. And for the 
strength of the material, it is very important whether the 
load sharing is global or local. Because, as we see from 
the phase diagram (see Fig. 0, the system breaks down 
as soon as any load is applied much before the load shar- 
ing is completely local. It is in fact close to (below) the 
percolation threshold in the respective lattices. 

Finally, we have compared our model with another in- 
terpolation scheme, which is in fact the homogeneous 
version of the present model. Although the average load 
sharing of local and global fractions are same, due to local 
fluctuations in the present model, the phase boundaries 
are different (see Fig. [T3"]) . 

In summary, we have proposed a heterogeneous load 
sharing mechanism for fiber bundle models to account 
for the local disorder concentration in the bottom plate 
of the fiber bundle models. For finite fraction of local 
load sharing, the model has loading history dependent re- 
sponse. We have found numerically, that the critical be- 
havior crosses over from GLS (mean field) to LLS at some 
effective site percolation threshold. The critical behavior 
across the phase boundary (Fig. [T3J) is of mean field type 
(relaxation time divergence exponent z = 0.50 ±0.01 and 
avalanche size distribution has an exponent 2.50 ± 0.01). 
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